Comment on "Superfluid stability in the BEC-BCS crossover" 



O 

o 

(N 



Daniel E. Sheehy and Leo Radzihovsky 
Department of Physics, University of Colorado, Boulder, CO, 
(Dated: August 7, 2006) 



80309 



We point out an error in recent work by Pao, Wu, and Yip [Phys. Rev. B 73, 132506 (2006)], that 
stems from their use of a necessary but not sufficient condition [positive compressibility (magnetic 
susceptibility) and superfluid stiffness] for the stability of the ground state of a polarized Fermi gas. 
As a result, for a range of detunings their proposed ground-state solution is a local maximum rather 
than a minimum of the ground state energy, which thereby invalidates their proposed phase diagram 
for resonantly interacting fermions under an imposed population difference. 
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There has been considerable recent interest in paired 
superfluidity of fermionic atomic gases under an imposed 
spin polarization^, i.e., when the numbers iV-j- and N± 
of the two atomic species undergoing pairing are differ- 
ent. Along with the detuning S of the Feshbach reso- 
nance (controlling the strength of the interatomic attrac- 
tion), the population difference AN = iVf — Ni is an 
experimentally-adjustable "knob" that allows the study 
of novel regimes of strongly-interacting fermions. 

A crucial question concerns the phase diagram of res- 
onantly interacting fermions as a function of S and im- 
posed AN. Two early theoretical studies that have ad- 
dressed this issue are the work by Pao, Wu, and Yip£ on 
the one-channel model of interacting fermions and our 
work^ on the two-channel model of interacting fermions, 
with the details and extensions of the latter presented 
in our recent preprint £ Apart from the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) phase appearing over a thin 
sliver of the phase diagram in Ref. |4| (explicitly not con- 
sidered by Pao, Wu, and Yip), for the positive-detuning 
BCS and crossover regimes the phase diagrams in these 
manuscripts are in qualitative agreement 6 . However, at 
negative detuning in the BEC regime, the phase dia- 
grams are qualitatively different. In particular, in Fig. 4 
of Ref. |3|, stable superfluidity is claimed to exist above 
and to the left of a nearly vertical and positively sloped 
phase boundary that we re- plot here in Fig. [T] as a dashed 
line at negative detuning. In qualitative contrast, in our 
work, Ref. 0, we found that stable superfluidity exists 
only below and to the left of a negatively-sloped phase 
boundary (see Fig. Q]) . 

What is the source of this qualitative discrepancy? 
Although Ref. and Ref. Q] use different models of 
resonantly-interacting fermions, the close relationship be- 
tween the one- and two-channel models (particularly 
within the mean-field approximation) implies that they 
should yield qualitatively similar phase diagrams. Fur- 
thermore, we have extended our original two-channel 
model study^ to that of a one-channel model^ and, as ex- 
pected, found results in qualitative agreement with those 
in our Letter—, but in disagreement with that of Ref. y- 

Indeed, as we explicitly show here, the origin of the 
discrepancy is that the criterion for stability of the super- 
fluid phase used in Ref.[3| (based on positivity of magnetic 
susceptibilities) is a necessary but not sufficient condi- 
tion for stability and does not ensurei^ that the state 
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FIG. 1: (Color online) Polarization AN/N vs. detuning 
— — — phase diagram of the one-channel model (appropriate 
for width 7 > 1) showing regimes of FFLO, superfluid (con- 
fined to the AiV = axis), magnetized superfluid (SFm ), and 
phase separation. The dashed lines are the purported phase 
boundaries reported in Ref. [jj. As we show^, for the entire 
shaded region single-component (uniform) solutions to the 
stationarity conditions Eqs. (J7J) are local maxima of Eg {A) 
(and therefore unstable, leading to phase separation), shown 
for the point marked with the "x" in Fig. [2] 



is even a local minimum of the energy. This thereby 
leads to incorrect phase boundaries in both the BEC and 
BCS regimes, although below we shall focus on the phase 
boundary inside the BEC regime, where the error is qual- 
itative and most pronounced. 

Thus, much of what is claimed to be a "stable super- 
fluid" in Ref. H (in the negative-detuning BEC regime) 
is actually unstable to phase separation. This is illus- 
trated in the correct T = mean-field phase diagram, 
Fig. [U for the one-channel model (quantitatively consis- 
tent with other recent wor k 10 ! 11 and derived in detail in 
Ref. [5|), plotted as a function of the dimensionless pa- 
rameter — (fepfls) 1 (proportional to the Feshbach res- 
onance detuning 5, with k-p the Fermi wavevector and 
a s the s-wave scattering length) and the polarization 
AN/N=(N r -N l )/(N 1 + N l ). 

Our starting point is the single-channel model Hamil- 
tonian (studied in Ref. [H) for two resonantly-interacting 
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species of fermion c k(T (a =f , J,): 



-k+q| L p-qT' 
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where e k = k 2 /2m, m is the fermion mass, and is 
the system volume. A broad Feshbach resonance is mod- 
eled by an attractive interaction A < 0, the magnitude 
of which increases with decreasing Feshbach resonance 
detuning. 

The equilibrium ground state of a many particle sys- 
tem in the grand-canonical ensemble (with chemical po- 
tentials jitf and /ii) at T = is characterized by the grand 
thermodynamic potential /zi ) defined byi£ 
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with the grand-canonical Hamiltonian 



H = H - MT^T - »l N l- 
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is the number operator for 



Here, N a 

fermion species a and the minimization in Eq. ([2]) is over 
all possible ground states. 

The standard (BCS-type) mean-field approximation 
that we shall utilize (as done in Ref. d) amounts to as- 
suming a restricted class of possible many-body ground- 
states self-consistently parametrized by the ground-state 
expectation value 



A 



(4) 



For pairing amplitude A ^ 0, this class includes both 
weakly-paired BCS-type and strongly-paired molecular 
BEC-type pairing order. For A = it is an unpaired 
Fermi gas. 

Once this mean-field approximation has been made, it 
is straightforward 4 * 5 - to compute the variational ground- 
state energy, which is the expectation value Eg(A) = 
(H). Evaluating this expectation value, and converting 
momentum sums to integrals, we find 



E G (A)/V = - 
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■{E k -h)G(h-E k ), 
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where = e k — fi, E k = \J £j[ + A 2 , and we have defined 
the chemical potential fi — ^(/if + /xjj and the chemi- 
cal potential difference h = h(fJ>i — Mi)- We have also 
exchanged the bare interaction parameter A for the vac- 
uum s-wave scattering length a s given by 
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The determination of the mean-field phase diagram us- 
ing Eq. ([5]) is conceptually quite simple. According to 



Eq. ([2]), the ground state at a particular fx-f and /ij (or, 
equivalently, /i and h) is given by the minimization of 
Eg (A) with respect to the pairing amplitude A that can 
be taken to be real. Any such minima of course satisfy 
the stationarity constraint or gap equation [equivalent to 
Eq. Q] 



= 



dE G 
dA 
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where we emphasize that the derivative is taken at fixed 
fi and h. Since experiments are conducted at fixed atom 
number, we must augment Eq. (|7a[) with the number con- 
straint equations N„ = (N a ) . Examining Eq. §5§ , we see 
that the constraints can be rewritten as 



N = 
AN = 



OEg 
dfi 
dE G 
' dh : 



(7b) 
(7c) 



with the total particle number N = + N± and popu- 
lation difference AN = iV T - N±. 

Our key point (apparently missed by the authors of 
Ref. 0) is that not every simultaneous solution of the gap 
and number equations, Eqs. (O, corresponds to a physi- 
cal ground state of the system; the additional criterion is 
that the solution A must also be a minimum of E G (A) 
at fixed The verification that an extremum solution 
is indeed a minimum is particularly essential when there 
is the possibility of a first-order transition, with Eg (A) 
exhibiting a local maximum that separates local minima, 
as is the case for a polarized Fermi gas, studied here and 
in Refs. iiH 

Analyzing Eqs. ([7]), we find that for sufficiently large 
AN in the positive-detuning BCS regime and for AN = 
N in the BEC regime, a solution to Eqs. ([7|) may be 
found that minimizes Eg (A) at A = 0, indicating a po- 
larized normal phase. Also, for sufficiently low detun- 
ing in the BEC regime, a polarized molecular superfluid 
(SFm ) solution exists that minimizes Eg(A) at A ^ 0. 
A more general analysi s 4 ' 5 ' 13 shows that a periodically- 
paired FFLO solution is the ground state over a thin 
range of polarization values at sufficiently large positive 
detuning. 

However, over the large shaded portion of the phase 
diagram, at intermediate detuning and polarization, we 
find that it is not possible to satisfy Eqs. ([7]) with a 
(homogeneous, single component) minimum of Eg(A). 
For the corresponding range of parameters the system 
phase-separates^ 4 - into two coexisting ground states (that 
are degenerate minima of Eg[A]). The resulting phase- 
separated state can be explicitly accounted for by gener- 
alizing the ground-state ansatz to include the possibility 
of such an inhomogeneous mixture 5 -. 

The contrasting strategy of Pao, Wu, and Yip 3 is to 
find solutions of Eqs. (O for all values of N, AN, and 
— k ^ a , some of which do not correspond to ground states. 
The unphysical (unstable) solutions are then discarded 
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FIG. 2: Plot of E G (A) at fixed fx = -0.013e F and h = 
1.16eF, normalized to the Fermi energy ef, for coupling 
g — fcga = 0.63, so that Eqs. ([7| yield a stationary solution 
with AN/N = 0.2 and A = 1.02. As seen in Fig. [2] however, 
this solution is a local maximum (saddle point) of Eg(A). At 
this point (^i- = 0.63 and AN/N = 0.2) in the phase di- 
agram (illustrated by an "x" in FigfTJ the actual mean-field 
ground state is a phase-separated mixture of a superfluid and 
a normal state. 



(thereby determining the phase boundaries plotted in 
Ref.S) based on criteria of the positivity of magnetic sus- 
ceptibilities (atomic compressibilities) and the superfluid 
stiffness. However, as we now discuss, these stability cri- 
teria are necessary but not sufficient (i.e., not generally 
restrictive enough) to ensure that a solution to Eqs. 
is indeed a minimum of Eq(A). In contrast, it can be 
shown that the converse is true, namely that the mini- 
mization of Eq(A) ensures that the local susceptibilities 
are positive definite^ 

The stability criteria^ used by Pao, Wu, and Yip can 
be understood by examining Fig. 3 of Ref. H The solid 
lines in this figure correspond to solutions of Eqs. ([7]) at 
different values of the coupling g = fcp g . In particular 
they plot h/ t-p (with e F the Fermi energy, related to the 

density n = N/V by n = |ce| /2 with c = to 3 / 2 / '^/2ir 2 h 3 ) 
as a function of the polarization nd/n = AN/N (at fixed 
density), where nd — AN /V is the magnetization. At 
positive and intermediate detunings (the bottom curves 
of Fig. 3 of Pao, Wu, and Yip), they find solutions satisfy- 
ing M^- 1 < and correctly conclude that such solutions 

° on d I n J 

(having a negative magnetic susceptibility) are unstable. 
However, at sufficiently negative detuning (g > 0.5) Pao, 
Wu, and Yip. find solutions to Eqs. (O with a positive 
susceptibility -§^-\ n > 0, and based on their criterion (er- 
roneously) conclude that these solutions indicate a sta- 
ble magnetic superfluid ground state. They then define a 
phase boundary in the BEC regime (the leftmost dashed 
curve of Fig. [J), to a stable magnetized superfluid, by 
where -M±- 1 changes sign. 

on d In ° ° 

However, our explicit calculation of -Eg (A) (plotted in 
Fig. (2) for one such solution (with = g = 0.63 and 
n d/ n — 0.2, indicated with an "x" in Fig. Q] and corre- 
sponding to a point on the uppermost solid curve of Fig. 3 



of Ref. d), purported by Pao, Wu, and Yip to be stable 
(to the left of their proposed stability boundary) , shows 
that in fact this solution (indicated with a dot in Fig. [2|) 
is a local maximum and therefore does not represent a 
ground state. This solution was obtained by numerically 
solving Eqs. ^ at g = 0.63 and AN/N = 0.2, yielding 
/i = -0.013e F , h = 1.16e F and A = 1.02e F , the latter 
two values consistent with Figs. 2 and 3 of Pao, Wu, 
and Yip, showing that we are indeed reproducing a so- 
lution claimed to be stable by Pao, Wu, and Yip. Thus, 
the method used by Pao, Wu, and Yip has not correctly 
located the global minimum of the ground-state energy; 
indeed, it has not even found a local minimum. 

Although it might appear from the plot of Eq(A, /i, h) 
(Fig. [2]) that the true ground state is an unpaired (A = 0) 
normal state, this state does not satisfy Eqs. (|7b[) and 
((7c)) ; thus, it is also not the ground state at this cou- 
pling and polarization. Indeed, as noted above, we find 
that it is impossible to minimize Eq(A) while satisfy- 
ing Eqs. ([7b]) and (fTcj) at this coupling and polarization 
(marked by an "x" in Fig. [T]), nor anywhere inside the 
shaded region in the phase diagram Fig. [TJ indicating 
the absence of a uniform solution. The true mean-field 
ground state everywhere in the shaded region is a phase- 
separated mixture of two phases of different densities in 
chemical equilibrium such that the total number and po- 
larization constraints are satisfied^. 

We note that a ground state determined by minimizing 
Eq{A) at a particular /if and /Xj automatically satisfies 
the condition of having a positive magnetic susceptibility 
(compressibility). Indeed, it is straightforward to gen- 
erally sho w 15 ' 17 that, since is concave down- 
wards^, the eigenvalues of the susceptibility matrix 



are positive in the ground state. An equivalent stability 
criterion was derived in Ref. [TH by considering stability 
against local density variations. The procedure used by 
Pao, Wu, and Yip, however, did not amount to analyzing 
Eq. ([H]); as discussed above, the phase diagram of Pao, 
Wu, and Yip was obtained by computing which is 

not an equivalent condition. If the authors "nad instead 
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studied 
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dh 



they would have found that the solution 
plotted in Fig.[2]has a negative magnetic susceptibility in 
the grand-canonical ensemble (and therefore is unstable). 

However, we must emphasize the important point that 
any particular extremum solution may have a positive 
magnetic susceptibility and still not be the ground state. 
The simplest example of this is the normal Fermi gas 
state (A = 0), which satisfies the gap and number- 
constraint equations everywhere in the phase diagram 
(including AA^ = 0) and has a positive magnetic suscep- 
tibility, but is only the actual ground state (a minimum 
of Eg{A)) at sufficiently large AA^. Thus, if the authors 
of Ref. |3j had computed the eigenvalues of Eq. ([5]) instead 
they would have been able to discard some of 
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the erroneous solutions plotted in Fig. 3 of Ref. 3. How- 
ever, in general, Eq. ||5J) is still not sufficient and the 
most correct scheme is to use Eq. ([2]), i.e., to find the 
global minimum, in the grand-canonical ensemble, of the 
mean-field ground-state energy. 

We conclude by noting that, although a mean-field 
analysis of the one-channel model is not expected to 
be quantitatively accurate near the resonance position 
where fcp|a s | — > oo, it is expected to yield a qualitatively 
correct description of a polarized resonantly-interacting 
Fermi gas. To summarize, we have shown that while for 



equal species number (AN = h = 0) such analysis can 
simply proceed by solving the gap and number equations 
[Eqs. 0], the existence of first-order transitions at h ^ 
implies that the ground-state energy Eq(A) exhibits local 
maxima as a function of A, yielding solutions to Eq. (|7a|) 
that do not represent physical ground states. 
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The second term of this equation is clearly positive at the 
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